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1 . 


STATEMENT OF THE PROBLEM 


a. Introduction 


The development of severe storms is influenced by inter- 
actions between scales of motion that range from planetary 
systems down to the microscale. Recent studies have focused on 
mechanisms that produce areas of intense convection, charac- 
teristics of the severe storm environment, and possible "feed- 
back" mechanisms from the convection to the larger-scale sur- 
roundings. Although feedback mechanisms are thought to modify 
synoptic-scale storm environments, they are not yet fully 
understood. Diagnostic studies of environmental kinematic 
parameters and the environmental kinetic energy balance are 
ways of examining scale interactions and feedback mechanisms in 
the severe storm environment. This report presents results 
from such a study. 

b. Previous studies 

1) Storm-environment interactions 

Interactions between convection and the surrounding larger- 
scale environment have been investigated for many years; 
however, most advances have been made only recently. Aubert 
(1957) used numerical models to show that latent heat release 
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associated with thunderstorms produced increased large-scale, 
low-level convergence and enhanced upper- level divergence. The 
synoptic-scale storm environment experienced low-level height 
falls and upper-level height rises. The effects of latent heat 
release on changes in pressure patterns have been described by 
Manabe (1956) . Latent heat release was found to enhance the 
production of kinetic energy (Danard, 1964, 1966) because 
induced ageostrophic flow was directed toward decreasing pres- 
sure in both upper and lower levels. Ninomiya (1971a and b) 
used conventional rawinsonde data and satellite imagery to 
detect a raid-tropospheric warm core, strong low-level conver- 
gence, and strong upper-level divergence near tornado producing 
thunderstorms. The formation of a mid-tropospheric jet stream 
near the storm area was attributed to the increased horizontal 
temperature gradient caused by the warm core. 

A recent surge of interest has developed concerning pos- 
sible feedback mechanisms from areas of large, and relatively 
long-lasting, intense storms. Based on enhanced infrared 
satellite imagery, Maddox (1980a and b) defined the Mesoscale 
Convective Complex (MCC) to be a nearly circular cloud shield 
encompassing intense convection lasting at least 6 h. The area 
of the shield exceeds that of typical individual thunderstorms 
by two orders of magnitude. Convectively driven, mesoscale 
circulations were hypothesized to be dominant within the MCC in 
contrast to larger scale features which control prefrontal 
squall lines. Environmental conditions associated with the 
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storm complexes have been described by Maddox (1979), Fritsch 
and Maddox (1980), Maddox et al. (1981), and Fritsch and Maddox 
(1981a). The most significant feature of the MCC was an area of 
mid-tropospheric warming and mean mesoscale ascent. In addi- 
tion, meso a-scale regions of high pressure and associated 
anticyclonic flow were observed in the upper troposphere. Jet 
maxima, whose winds often were 20 m s ^ greater than those pre- 
dicted by the National Weather Service’s Limited Fine Mesh Model 
(LFM) , formed on the poleward side of the complexes. Synoptic- 
scale divergence over the storm area was typically 15 x 10 ^ 
s ^ . These environmental changes were attributed to the storms; 
however, the exact physical mechanisms involved require further 
study. A significant point of the investigations was that 
current numerical forecasting models do not effectively simulate 
the major large-scale modifications associated with MCC's be- 
cause of our incomplete understanding of their complex nature. 
Other recent diagnostic studies that have described mean, meso- 
scale ascent in response to convective storms include those by 
Fankhauser (1969, 1974), Sanders and Paine (1975), Sanders and 
Emanuel (1977), and Ogura and Chen (1977). It is noteworthy 
that numerical simulations have duplicated some of the modifi- 
cations attributed to severe storms (e.g., Rao and Hassebrock, 
1972; Kreitzberg and Perkey, 1977; Fritsch and Chapell, 1980a 
and b; Fritsch and Maddox, 1981b); however, a great deal more 
remains to be learned before the effects of large, intense areas 
of convection can be included in the prediction models. 
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2) Kinetic energy analyses of the total flow 


Energy levels of the general circulation, together with the 
various sources, sinks, and transformations of energy, have been 
studied extensively (e.g., Dutton and Johnson, 1967; Newell et 
al. , 1970; Saltzman, 1970; Peixbto and Oort, 1974; Oort and 
Peix6to, 1974). Other investigations have examined the ener- 
getics of more limited regions, especially extratropical cy- 
clones. Smith (1980) presents an excellent review of cyclone 
energetics in which he notes that such systems are energetically 
active components of the general circulation. Individual cy- 
clones often were found to exhibit major departures from the 
average . 

Large areas of intense convection can be associated with 
cyclone activity or, in the case of MCC's, can occur without 
major large-scale forcing mechanisms. In either case, the 
larger-scale energetics of the convective environment have 
received relatively little attention. Fuelberg and Scoggins 
(1978) described the synoptic-scale kinetic energy balance of 
two large storm areas occurring during NASA's fourth Atmospheric 
Variability Experiment (AVE 4, 24-25 April 1975). Although no 
major cyclone activity was present, the various kinetic energy 
transformations and transports were quite intense. In fact, the 
storm areas appeared to influence the general circulation as 
much or more per unit area than intense cyclones. With the aid 
of the special 3 h rawinsonde data, they observed major fluctua- 
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tions in energy balance that seemed related to the life cycles 
of the intense convection within the region. Similarly, 
Robertson and Smith (1980) found major kinetic energy activity 
in the extratropical cyclones that produced the Jumbo (3-4 April 
1974) and Palm Sunday (10-11 April 1965) tornado outbreaks. 
Companion studies determined that latent heat release associated 
with intense convection provided the greatest source of avail- 
able potential energy to the areas (Lin and Smith, 1979; Smith 
and Lin, 1980). An intriguing new hypothesis is that widespread 
convection can act as a source of synoptic-scale kinetic energy. 
Studies by Kornegay and Vincent (1976), Sheu and Smith (1977), 
and Vincent and Schlatter (1979) suggest that convection plays a 
role in the transfer of kinetic energy from unresolvable (sub- 
grid) to resolvable (grid) scales of motion which appears in the 
"dissipation" term of the kinetic energy balance equation. 

Very few kinetic energy studies have been conducted on the 
subsynoptic scale, mainly due to the lack of appropriate data. 
Mclnnis and Kung (1972), Kung and Tsui (1975), and Tsui and Kung 
(1977) used meso 3 -scale rawinsonde data (station spacing of 80 
km) from the National Severe Storms Laboratory (NSSL) in their 
investigations. Some energy transformations on this scale were 
greater than those observed in major synoptic-scale cyclones. 

The magnitudes of mesoscale energy processes were found to vary 
greatly, depending on the strength of the convective systems in 
the area. Time variations in the energy budget terms seemed 
related to the growth and decay of the nearby storms. 
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Data from the recent SESAME periods will allow descriptions 
of mesoscale atmospheric phenomena since the average spacing of 
the rawinsonde network was either 250 km (regional scale) or 100 
km (storm scale). Carlson et al. (1980) and Moore and Fuelberg 
(1981) have described the rapid changes in horizontal wind, 
divergence, vertical motion, and stability that occurred in 
conjunction with the deadly Red River Valley tornado outbreak of 
10-11 April 1979. Fuelberg et al. (1980) described the kinetic 
energy balance of the synoptic-scale storm environment based on 
National Weather Service (NWS) sonde data alone, while Jedlovec 
and Fuelberg (1981) performed a similar analysis using both NWS 
and special site data such that meso a-scale resolution was 
obtained. Area averaged energetics at individual observation 
times showed that the energy balance during periods of maximum 
convective activity differed considerably from that of the 
remaining periods. The local kinetic energy balance over 
Oklahoma during the formation of a limited upper-level wind 
maximum received special attention. Cross-contour production of 
kinetic energy was found to be the dominant local source for its 
development. This raised the possibility that the intense 
convection associated with the Red River Valley Outbreak may 
have contributed to the formation of the wind maximum. 
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3) Kinetic energy analyses of divergent and rotational 
components 

A different approach to studying storm-environment inter- 
actions is to evaluate the energetics of the divergent and 
rotational wind components, and V^, respectively. Such an 
approach may further explain the mechanisms by which intense 
thunderstorms modify their larger-scale environments and how 
such alterations might be included in the numerical forecast 
models. Diagnostic studies of this type are quite rare, 
however. 

Chen and Wiin-Nielson (1976) have shown that available 
potential energy is converted into kinetic energy of the diver- 
gent wind (Kp) and then into kinetic energy of the rotational 
component (K^) . Quantity was found to be very important 
because all of the converted available potential energy goes 
through Kp which remains small and thereby plays a catalytic 
role in the conversion process. Northern Hemispheric summertime 
data showed that K^/Kj. = 10% (Kp, is the total kinetic energy) 
with a maximum near the level of the jet stream in middle lati- 
tudes. Krishnamurti (1971) and Chen (1980) showed that is 
the dominant part of the flow at 200 mb in the tropics. Chen et 
al. (1978), who studied a major mid-latitude cyclone, found that 
although Kpj/Kp = 2.5%, assumed crucial importance in the 
generation and horizontal transport of kinetic energy. Simi- 
larly, Smith (1974) noted that alterations to the divergent wind 
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produced significant changes in station values of total kinetic 
energy generation and horizontal flux divergence. On the other 
hand, Krishnamurti (1968) determined that much of the cross- 
contour flow near a developing wave cyclone could be explained 
by V R and was attributable to differential vorticity advection 
and the Laplacian of thermal advection. 

c . Objectives 

The overall goal of this research is to investigate the 
roles of the divergent and rotational wind components in the 
kinetic energy budget of a severe storm environment. A diag- 
nostic study is performed to obtain the synoptic-scale kinetic 
energy balance of two large storm areas that occurred during 
NASA's AVE 4 period (24-25 April 1975). Synoptic conditions 
during AVE 4 and the energetics of the total flow have been 
studied previously (e.g., Fuelberg and Scoggins, 1978; Maddox, 
1980 a and b; Maddox et al. , 1981). Apparently, however, rela- 
tive contributions of the divergent and rotational wind com- 
ponents to total kinetic energy have not been evaluated during 
any period of intense convection. The following points will be 
investigated: 

1) Space variabilities of V R , V^, K R and K^, 

2) Time variabilities of these quantities, 
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Contributions of and to the kinetic energy 
balance, and 

Error bounds for all derived results. 


2 . 


THEORETICAL CONSIDERATIONS 


The kinetic energy equation for a fixed, limited volume in 

isobaric coordinates is given by Smith (1969) as: 

0K_/3t = ff -V* Vcf> - ff V-kT/ - ff Scok/Sp + D + / k 9p /8t , (1) 

o o 

(A) (B) (C) (D) (E) 

where 

ff = 1/g fff dx dy dp, 

-> 

V is the horizontal wind vector, 
w is the vertical motion in isobaric coordinates, 
k = (u + v )/2 is the horizontal kinetic energy per 
unit mass, 

Kp = ff k, 

<j> = gz is the geopotential height, 

A is the computational area, and 
subscript o denotes surface values. 

On the right hand side of (1), term A represents generation 
of kinetic energy due to cross-contour flow while terms B and C 
represent horizontal and vertical flux divergence of kinetic 
energy, respectively. Term D is computed as a residual to (1) 
and represents both thermodynamic and mechanical frictional 
processes as well as transfers of energy between resolvable and 
unresolvable scales of motion. This term also contains possible 
errors arising from all other terms of the equation and is 
commonly referred to as the "dissipation” term. Term E repre- 
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sents variations in kinetic energy due to changes in the mass of 
the volume under consideration. Since it is several orders of 
magnitude smaller than the other terms in (1) , it will not be 
considered further. In regions of large pressure change, the 
term could be significant; however, this is not observed during 
the AVE 4 period. 

To examine contributions of the divergent and rotational 
wind components to the total kinetic energy balance, these 
components must be introduced into (1). A summary of the Chen 
et al. (1978) derivation for that purpose follows. 

.-V. 

Helmholtz’s theorem states that the horizontal velocity (V) 
can be expressed as the sum of the divergent wind (Vjj) and the 
nondivergent or rotational wind 0^) , i.e., 

* = + V < 2) 

Kinetic energy per unit mass is given by 

k = % f = kp + k D + tf R > ? D , (3) 

where 

^ % V * R and kjj = * • V (4) 

Using (2) , the integrands of terms A and B in (1) can 
be written as 

, (5) 

. -V-k^ R - ^'kV D , (6) 

Equations (5) and (6) show that horizontal flux divergence and 
generation of kinetic energy, respectively, are due to both the 
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divergent and nondivergent wind components . Terms -V * V<}> and 

K. 

_V ' often are denoted as the generation of kinetic energy 

due to barotropic and baroclinic processes, respectively (Chen 
and Wiin-Nielson, 1976). Using (5) and (6), (1) can be re- 
written as 

SK^/at = ff -V R - ^ + ff -V D * V<f> - ff ^*kV R 

- ff ^‘kV D - ff 9tok/3p + D . (7) 

The emphasis of this research is to describe the relative contri- 
butions of and to kinetic energy content in (3) , horizontal 
flux divergence in (5), and generation of kinetic energy in (6). 
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3. DATA AND COMPUTATIONS 


a. Data 


The case being investigated is the National Aeronautics and 
Space Administration's (NASA) fourth Atmospheric Variability 
Experiment (AVE 4), conducted on 24-25 April 1975. The AVE 4 
period has been widely studied because it contains two severe 
storm outbreaks (e.g., Wilson, 1976; Fuelberg and Scoggins, 

1978; Maddox, 1980a and b; Maddox et al., 1981). Rawinsonde 
data from 42 National Weather Service (NWS) stations were taken 
at nine times - 24 April at 0000 GMT, 0600 GMT, 1200 GMT, 1500 
GMT, 1800 GMT, and 2100 GMT, and on 25 April at 0000 GMT, 0600 
GMT, and 1200 GMT. Locations of these stations are shown in 
Fig. 1. The rawinsonde data were given at 25 mb intervals by 
Fucik and Turner (1975) while the data reduction procedures used 
on this set have been described by Fuelberg (1974). 

All rawinsonde soundings were checked carefully to deter- 
mine if sondes had entered violent thunderstorm updrafts or 
downdrafts. Nonhydrostatic accelerations such as thunderstorm 
updrafts are a potential source of error since hydrostatic 
equilibrium is assumed. Data for four soundings (HTS - 1800 
GMT; TOP - 0600 GMT; UMN - 1200 GMT; and DAY - 0600 GMT, 24 
April) that experienced severe drafts were removed from the data 
set. 
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Fig. 1. 


Ravinsonde stations 
AVE 4 experiment. 


participating in the 
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Surface data from 310 stations were obtained from the 
National Climatic Center and used in the investigation. 

Manually Digitized Radar (MDR) data were obtained from the 
National Oceanic and Atmospheric Administration's (NOAA) Tech- 
niques Development Laboratory to determine objectively the 
intensity and location of the radar-observed convection. During 
1975 MDR data were coded with a single digit from 0 to 9 to 
indicate areal coverage and echo intensity within blocks that 
were approximately 83 km on a side. Table 1 defines the MDR 
code. Plots of MDR values were made each hour for 3 h periods 
centered on each of the nine rawinsonde observation times. The 
three hourly plots then were combined into a single chart for 
each of the nine times by using the highest coded value reported 
for individual MDR blocks. 

b . Analytical procedures 

1) Objective analysis 

An objective analysis scheme by Barnes (1964) was used to 
interpolate data from the randomly spaced rawinsonde and surface 
stations to equally spaced grid points. The grid network, 
centered over the AVE 4 area, had a spacing of 158 km. Data 
from each rawinsonde site influenced grid points within a scan 
radius of three grid distances while four iterations were allowed. 
The surface data were analyzed using a scan radius of two grid 
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Table 1. Manually digitized radar (MDR) data code 
(Foster and Reap, 1973). 


Code No. 

Maximum 
Observed 
VIP Values 

Coverage 
In Box 

Maximum 
Rainfall, 
Rate (in h i ) 

Intensity 

Category 

0 

No Echoes 




1 

1 

Any VIP1 

< 0.1 

Weak 

2 

2 

£ 50% of VIP2 

0. 1-0.5 

Moderate 

3 

2 

> 50% of VIP2 

0. 5-1.0 

Moderate 

4 

3 

£ 50% of VIP3 

1. 0-2.0 

Strong 

5 

3 

> 50% of VIP3 

1. 0-2.0 

Strong 

6 

4 

£ 50% of VIP3 
and 4 

1. 0-2.0 

Very Strong 

7 

4 

> 50% of VIP3 
and 4 

1. 0-2.0 

Very Strong 

8 

5 or 6 

£50% of VIP3, 
4, 5, and 6 

> 2.0 

Intense or 
Extreme 

9 

5 or 6 

> 50% of VIP3, 
4, 5, and 6 

> 2.0 

Intense or 
Extreme 


distances with four iterations. Even though the resulting 
objectively analyzed fields agreed well with subjectively (hand- 
analyzed) data fields, a mild filter by Shuman (1957) was applied 
to the grids to further suppress features below the minimum data 
resolution (twice the station separation) . The smoothing element 
index used in this procedure was 0.2. 

Gridded analyses of height, temperature, wind components, 
and moisture content were produced at 18 levels, i.e., the 
surface and at 50 mb intervals from 900 mb to 100 mb for each of 
the nine times. Winds at the 18 levels were arithmetic averages 
over 50 mb layers to reduce the effects of random errors that 
are inherent in the data. The grids resulting from the above 
procedures are identical to those employed by Wilson (1976) and 
Fuelberg and Scoggins (1978). 

2) Computation of divergent and nondivergent winds 

An iterative scheme developed by Endlich (1967) and used by 
Chen et al. (1978) was employed to separate the horizontal wind 
into its divergent and rotational components . This scheme 
allows the divergence to be virtually eliminated while con- 
serving the original vorticity. A brief description of this 
method follows. 

First, vorticity is calculated using centered finite dif- 
ferences on the original analyzed wind components. To reduce 
the divergence at a given grid point to zero, adjustments of 
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equal magnitude, but opposite sign, are made to the u and v wind 
components at the neighboring grid points. While reducing the 
divergence, the procedure modifies the vorticity of the resulting 
wind. Therefore, to retain the original vorticity, u and v in 
the vicinity of the grid point in question are modified appro- 
priately. Then at the next grid point, the most recent winds 
are used. The procedure continues until all grid points are 
handled; then the second iteration begins. With repeated appli- 
cation of this procedure, the divergence of the adjusted wind 
approaches zero while its vorticity remains that of the original 
wind. For our purposes, the divergence of the adjusted wind 
field was everywhere reduced to at least 1 x 10 ^ s which is 
less than 0.1% of the largest magnitudes in the original winds 
(1.4 x 10 ^ s ^). Thus, the final field can be considered to be 
nondivergent in character. A constant correction was added to 
the u and v components of the nondivergent winds so that average 
values agreed with those of the original components. Similar 
corrections were made by Endlich (1967) and Chen et al. (1978). 
Finally, the divergent wind was obtained as the difference 
between components of the nondivergent and the original winds. 

An alternative technique for computing rotational and di- 
vergent winds is the solution of Poisson equations for stream 
function and velocity potential (Charney, 1962). This procedure 
was investigated for use here. The need for specifying meteoro- 
logically satisfactory boundary conditions in solving the Poisson 
equations meant a considerable loss of area on the periphery of 
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the AVE 4 region. Also, computational limitations initially 
resulted in a 30% loss of divergence when compared with values 
obtained from the original winds. Although later use of a 
staggered grid network was successful in retaining original 
divergence, it required that the divergent and nondivergent 
winds be interpolated back onto the original grid before solu- 
tion of the kinetic energy budget equation. Because of its 
simplicity and accurate results, the Endlich procedure was 
utilized in this study. 

3) Computation of vertical motion 

The kinematic method was used for computing vertical mo- 
tion. As the lower boundary condition, a value of zero was 
assumed at the surface. A limitation of the kinematic method is 
that errors in the computation of divergence accumulate by the 
end of the integration process, causing unrealistic values of 
vertical motion near 100 mb. To reduce this problem, the scheme 
suggested by O’Brien (1970) was used to adjust the profile of m 
to zero at 100 mb. Adjusted values of divergence were not used 
in the calculation of V^. 

4) Finite differencing and integration 

Centered finite differences were used to compute horizontal 
derivatives, and all vertical derivatives except those at the 
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surface and at 100 mb where noncentered differences were 
employed. Time derivatives also were evaluated using centered 
differences where possible, but uncentered differences at the 
first and last time were required. 

Terms of the kinetic energy budget equation (7), except 
dissipation, were computed at each grid point for each of the 18 
levels and then integrated over 50 mb layers from the surface to 
100 mb using the trapezoidal rule. The dissipation term was 
computed as a residual to balance (7) at each grid point in each 
50 mb layer. 

c. Error analysis 

Rawinsonde data contain errors which influence all sub- 
sequent calculations to varying degrees (Kurihara, 1961). To 
obtain quantitative confidence limits for the kinetic energy 
budget terms of the present study, separate computations were 
performed using randomly perturbed data. The procedures uti- 
lized in this analysis, as well as detailed results, are pre- 
sented in the Appendix. Generally, however, errors in rawin- 
sonde data should not seriously affect the interpretations of 
results that follow. One should exercise some discretion in 
assigning significance to minor term fluctuations, especially in 
the upper levels. 
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4. 


SYNOPTIC CONDITIONS 


Although the AVE 4 period was characterized by weak pres- 
sure and frontal activity, two large areas of intense convection 
occurred. Both of these areas have been classified as Mesoscale 
Convective Complexes (MCC’s) (Maddox, 1979). By the end of the 
36 h period, eight people had died and at least fifty injured as 
a result of a dozen tornadoes, wind damage, and flash flooding. 

Since detailed descriptions of weather conditions during 
AVE 4 have been given by Fuelberg and Scoggins (1978) , Maddox 
(1979), and Maddox et al. (1981), only a brief discussion of the 
major synoptic-scale features will be presented here. At 0000 
GMT 24 April (Fig. 2), the upper troposphere was characterized 
by a poorly defined polar jet stream over the northern United 
States and a subtropical jet along the Gulf Coast States. 

Intense thunderstorms (Fig. 3) were occurring over Kansas, 
Missouri, and Illinois in response to a weak short-wave trough 
in the middle troposphere (Fig. 2). Strongest thunderstorms 
were located at the nose of a weak speed maximum in the polar 
jet stream. No pronounced diffluence was evident over the AVE 4 
area at this time. 

By 0600 GMT (not shown), definite changes had occurred on 
the synoptic scale. The thunderstorms had organized into the 
first MCC of the period (Fig. 3). Rapid storm development had 
occurred eastward of the inital locations as the short-wave 
trough at 500 mb advanced. Additional formation occurred over 
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(b) 850 mb. 

Synoptic conditions at 0000 GMT 24 April 1975 
(Fucik and Turner, 1975). 
















National Meteorological Center radar summaries 
for the first storm outbreak. Echo tops are in 
kilometers . 






Kansas. At jet stream level, the polar jet became much stronger 

along the northern edge of the convection. Speeds had increased 
**1 

as much as 15 m s during the previous 6 h period over Illinois 
and Indiana, and now were as great as 65 m s ^ . Speeds had 
decreased south of the MCC. Pronounced diffluence in the polar 
jet now was evident over northern Missouri, Indiana, and Ohio. 
These significant changes in synoptic-scale flow patterns are 
hypothesized to be due to feedback processes from the MCC’s to 
the synoptic scale (Maddox, 1979). Thunderstorm activity re- 
mained intense at 1200 GMT (Fig. 3), but the storms weakened to 
showers by 1500 GMT (not shown) as they moved to the east. 

The second MCC initially developed near 2100 GMT. The 
outbreak began in eastern Kansas along a cold front extending 
from a low in Kansas, southwestward into Oklahoma (Fig. 4). At 
500 mb, a second short-wave trough over the Dakotas caused the 
surface cyclone in central Kansas to intensify. In the upper 
troposphere, the subtropical jet was located over Texas and 
Louisiana with a strong speed maximum over Texas. The polar jet 
curved from the Texas Panhandle into the northeastern United 
States. Strongest winds were over Nebraska, with secondary 
maxima over Illinois and southwest Missouri. A broad band of 
weakly defined diffluence covered the southeast half of the AVE 
4 area. By 0000 GMT 25 April (Fig. 5), a storm complex stretched 
from Missouri into Oklahoma with maximum radar tops reaching 
55,000 ft (16.8 km). 
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By 0600 GMT (Fig. 5) the newly formed MCC was located over 
the Mid-Mississippi Valley. Strongest thunderstorms stretched 
from eastern Kentucky into eastern Oklahoma, ahead of the 
short-wave trough at 500 mb. In the upper troposphere (not 
shown) , remarkable changes again had occurred in the flow 
patterns. Although peak speeds along the subtropical jet had 
decreased 10 to 20 m s \ a wind maximum again had developed on 
the polar jet just north of the strongest radar echoes. 

Greatest speeds within this area exceeded 65 m s \ an increase 
of over 15 m s 1 from 6 h earlier. The upper-tropospheric flow 
now was extremely diffluent over the region of strong thunder- 
storm activity; peak values of divergence at 200 mb reached 
14 x 10~ 5 s’ 1 . 

The MCC had moved southeastward by 1200 GMT (Fig. 5). 

Radar echoes extended along an arc from eastern Kentucky, into 
northern Alabama and Arkansas. The polar jet stream remained 
intense and continued to exhibit strongly diffluent flow over 
the active storm area (Fig. 6). Its location continued to be 
along the northern edge of the cirrus shield associated with 
the thunderstorm complex. 
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(d) 300 mb 


(Continued) . 





5. 


RESULTS 


The energy balance was computed for each observation time 
over the entire AVE 4 region (Fig. 1). Results for the single 
times then were averaged over the entire 36 h period. Energy 
budgets also were compiled for smaller, more limited volumes 
that moved with the convection. This section begins by giving 
an overview of the energetics for the entire region. Later 
sections focus on the energy balances of the smaller subregions. 
Since Fuelberg and Scoggins (1978) have already described the 
kinetic energy balance of the total flow, this discussion empha- 
sizes roles of the divergent and nondivergent wind components. 

a. Energetics of the composite AVE 4 period 

Average energetics of the large-scale flow in which the 
storms are imbedded is described by the composite budget for all 
nine observation times. Table 2 gives values for the total 
flow. Greatest energy content is found between 300-200 mb in 
association with jet stream activity, but total kinetic energy 
decreases with time at all levels during the 36 h period. 
Dissipation to subgrid-scale motions is a major factor producing 
this decrease below 800 mb where surface friction dominates. 
However, destruction of kinetic energy due to supergeos trophic 
cross-contour flow is the largest energy sink in the upper 
atmosphere (above 600 mb). Vertical flux divergence is found in 
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Table 2. Area-time average^kinetic energy budget. 

All units are W m except for IL, which is 
1(T Jm 


Pressure 
Layer (mb) 

S: 


-V T - Vif 

V*kV T 

3u)k/3p 

D 

200-100 

4.40 

-1.16 

-2.98 

0.05 

-1.42 

0.45 

300-200 

5.41 

-0.52 

-2.24 

1.80 

-1.47 

2.05 

400-300 

3.37 

-0.25 

-1.47 

-0.35 

0.27 

1.14 

500-400 

2.16 

-0.29 

-0.81 

-0.28 

0.29 

0.53 

600-500 

1.54 

-0.37 

-0.40 

0.09 

0.44 

0.56 

700-600 

1.07 

-0.37 

0.19 

0.29 

0.63 

0.36 

800-700 

0.79 

-0.33 

0.13 

0.09 

0.53 

0.16 

900-800 

0.70 

-0.44 

1.40 

-0.05 

0.48 

-1.41 

SFC-900 

0.35 

-0.24 

2.21 

0.01 

0.25 

-2.19 

Vertical 

Total 

19.8 

-4.0 

-4.0 

1.6 

0 

1.6 



the lower layers while flux convergence occurs above, thus 
indicating an upward transport of energy by rising air asso- 
ciated with widespread convection. Horizontal flux divergence 
of kinetic energy is small through most of the atmosphere except 
near 250 mb where it provides an important sink. (Here, the 
sign convention is that positive values of vertical and hori- 
zontal flux indicate an outward transport of energy). Positive 
dissipation is an important energy source above 800 mb, sug- 
gesting transfers from subgrid to grid scales of motion. 

Table 3 shows the distribution of kinetic energy between 

the divergent and rotational wind components (K^, K R ) and term 

V R . V Values of the ratio K^/K^, also are given. Maximum 

values of Kp, K R and K^ occur in the 300-200 mb layer, thereby, 

coinciding with the jet stream. Values of K R are only slightly 

smaller than those of Kp. The ratio K^/Kp is small above 500 mb 

with a maximum slightly greater than 2%. Below 500 mb, however, 

values of K^/Kp increase toward the surface as magnitudes of Kp 

and K r diminish. The error in approximating Kp by K R is always 

->• -y 

less than 10% during the composite period. Term V R * V^ is smaller 
in magnitude than K^ at all levels and is negative above 900 mb. 
Its sign and magnitude are dependent on the orientation of the 
two wind vectors, e.g., negative values occur when the winds 
oppose each other. Although the term integrates to zero in a 
global domain, this does not occur when limited regions such as 
AVE 4 are considered. 


34 



Table 3. Components of area-time averaged^, kinetjc energy 
content. Units of energy are 10 J J m 


Pressure 
Layer (mb) 


*R 

*D 

V • V 
R D 


200-100 

Oo 

4.36 

0 . 06 “ 

- 0.02 


300-200 

5.41 

5.35 

0.12 

-0.06 

2 . 2 % 

400-300 

3.37 

3.32 

0.05 

0.00 

1.5% 

500-400 

2.16 

2.15 

0.03 

- 0.02 

1 . 6 % 

600-500 

1.54 

1.53 

0.04 

-0.03 

2 . 6 % 

700-600 

1.07 

1.05 

0.04 

- 0.02 

3.7% 

800-700 

0.79 

0.78 

0.04 

-0.03 

5.0% 

900-800 

0.70 

0.66 

0.05 

- 0.01 

7.0% 

SFC-900 

0.35 

0.31 

0.03 

0.01 

8 . 6 % 

Vertical 

Total 

19.79 

19.51 

0.46 

-0.18 

2.3% 


Horizontal transport of kinetic energy into and out of the 
region is represented by the flux divergence term ( V*kV ). 

Table 4 presents a vertical tabulation of the divergent and 
rotational components comprising the total flux. The divergent 
wind is a source of flux convergence below 700 mb but contri- 
butes to flux divergence above 600 mb. Greatest values of 
V’kVp occur at jet stream level. For the most part, hori- 
zontal flux due to the rotational wind opposes the contribution 
of the divergent wind. Term V-kV R is a source of flux 
convergence above 500 mb and flux divergence nearer the surface. 
Both components yield flux divergence in the midtroposphere. 
Although Table 3 revealed that is one to two orders of magni- 
tude smaller than K R , the divergent wind component is as important 
as the rotational component in contributing to the total horizontal 
flux. In fact, the divergent component dominates V*kV^ above 
300 mb. The vertically integrated result (Table 4) is that the 
divergent component provides a major sink of energy to the 
region, but 64% of this amount is offset by inflow due to the 
rotational wind. 

Components of the cross-contour generation term are given 
in Table 4 as a function of pressure. As observed with the flux 
term, the divergent and rotational components produce differing 
results above 700 mb. The divergent component crosses the 
contours toward lower values at most levels, thereby generating 
kinetic energy. The rotational component produces the opposite 
result above 700 mb. The combined effect is destruction of 
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Table 4. 

Components of area-time averaged generation and 
horizontal flu^ divergence of kinetic energy. 
Units are W m . 



Pressure 
Layer (mb) 


-4- -> 

-V R *V# 

-%-h 

— y -V 

V*kV T 


-4- -4- 

V-kV D 

200-100 

-2.98 

-3.01 

0.03 

0.05 

-1.45 

1.50 

300-200 

-2.24 

-2.93 

0.69 

1.80 

-0.78 

2.58 

400-300 

-1.47 

-2.13 

0.66 

-0.35 

-0.87 

0.52 

500-400 

-0.81 

-1.11 

0.30 

-0.28 

-0.39 

0.11 

600-500 

-0.40 

-0.89 

0.49 

0.09 

0.05 

0.04 

700-600 

0.19 

-0.25 

0.44 

0.29 

0.29 

0.00 

800-700 

0.13 

0.26 

-0.13 

0.09 

0.20 

-0.11 

900-800 

1.40 

1.18 

0.22 

-0.05 

0.11 

-0.16 

SFC-900 

2.21 

1.21 

1.00 

0.01 

0.08 

-0.07 

Vertical 

Total 

-4.0 

-7.7 

3.7 

1.6 

-2.8 

4.4 


kinetic energy above 600 mb; however, in the lower troposphere 
(below 800 mb), both components produce kinetic energy. The 
vertical totals indicate that generation by offsets by 48% 
the destruction of kinetic energy by V^. It is significant that 
divergent generation provides the only source of kinetic energy 
to the AVE 4 region during the composite 36 h period. 

b. Time variability for the entire area 

The special AVE data set allows energy changes at 3 and 6 h 
intervals to be monitored during the convectively active period. 
Time variability of the surface-100 mb total flow energy balance 
is given in Table 5. Since boundary conditions for vertical 
velocity are zero at both the surface and 100 mb, vertical flux 
in the column integrates to zero, and therefore is not included. 
Total kinetic energy decreases with time except at 1200 GMT and 
2100 GMT 24 April when slight increases are noted. The time 
derivative reflects these variations; however, one should recall 
that centered time differences were used in the computations. 
Horizontal flux convergence is an energy source at the start of 
AVE 4 and from 1200 GMT to 2100 GMT 24 April. Strong flux 
divergence occurs near peak intensities of both MCC’s (0600 GMT 
24 April and 0600 GMT 25 April) as well as near the end of the 
experiment. 

Destruction of kinetic energy by cross-contour flow pro- 
vides a major sink between 0600 GMT and 1500 GMT 24 April which 
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Table 5. Vertical totals of the area averaged kinetic energy 
budget fgr individual observation timtjjs. All units 
are W m except Kp which is 10 J m 


Date/Time 

Kt 

BKj/at 

- v T * V < j > 

->* 

V'kV T 

D 

24/00 

23.01 

- 11.39 

0.77 

- 6.07 

- 18.23 

24/06 

20.55 

- 5.40 

- 10.48 

8.68 

13.76 

24/12 

20.70 

- 3.82 

- 24.26 

- 1.53 

18.91 

24/15 

19.31 

- 6.37 

- 12.77 

- 2.43 

3.97 

24/18 

19.31 

1.52 

9.87 

- 3.35 

- 11.70 

24/21 

19.65 

- 1.17 

5.61 

- 2.17 

- 8.95 

25/00 

19.07 

- 2.29 

- 5.84 

1.03 

4.58 

25/06 

18.86 

- 2.55 

2.72 

4.68 

- 0.59 

25/12 

18.00 

- 4.39 

- 1.50 

15.85 

12.96 

Composite 

19.8 

- 4.0 

4.0 

1.6 

1.6 



includes much of the life cycle of the first MCC. Destruction 
is due to supergeostrophic flow being directed across the con- 
tours toward higher pressure. Although weak generation occurs 
near peak intensity of the second MCC (0600 GMT 25 April) , 
strongest values are found between 1800 GMT and 2100 GMT. Signs 
of dissipation and generation are opposite at every observation 
time. During times of cross-contour destruction, kinetic energy 
is transferred from subgrid to resolvable scales (positive 
dissipation) and vice versa. Positive dissipation provides a 
major source of synoptic-scale energy between 0600 GMT and 1500 
GMT 24 April as well as at 1200 GMT 25 April. 

It is interesting to note that there is little correlation 
between the kinetic energy balance of the total flow and the 
life cycles of the two areas of intense convection, at least 
when the entire region is considered. When budgets for limited 
volumes enclosing the two MCC's are described in a later sec- 
tion, a clearer relationship will be apparent. 

Time variability of kinetic energy partitioned among the 
divergent and rotational wind components is given in Table 6 for 
the surface-100 mb layer. Rotational energy remains slightly 
smaller than K,j, until the last observation time, and temporal 
variations of closely reflect those of K^,. On the other 
hand, there are two maxima in the divergent component, and both 
occur at times of peak convective activity. Values of 
reveal that comprises between 1.6% and 2.7% of the total 
energy during periods of little convection. During times of 
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Table 6. Vertical totals of the area averaged components of 
kinetic energy content for ind^viduaj observation 
times. Units of energy are 10 J m 


Date /Time 



*D 

-* -> 

V -V 
R D 


24/00 

23.01 

22.54 

0.36 

0.11 

or 

24/06 

20.55 

19.76 

0.70 

0.09 

3.4% 

24/12 

20.70 

20.08 

0.48 

0.14 

2.3% 

24/15 

19.31 

19.03 

0.39 

-0.11 

2.0% 

24/18 

19.31 

19.21 

0.35 

-0.25 

1.8% 

24/21 

19.65 

19.63 

0.37 

-0.35 

1.9% 

25/00 

19.07 

18.98 

0.33 

-0.24 

1.7% 

25/06 

18.86 

18.30 

0.85 

-0.29 

4.5% 

25/12 

18.00 

18.03 

0.48 

-0.51 

2.7% 

Composite 

19.8 

19.5 

0.5 

-0.2 

2.6% 



peak MCC activity, however, values of 3.4% and 4.5% are observed 
for the first and second outbreaks, respectively. Thus, the 
enhanced divergence associated with the convective storms is 
evident even when the entire region is considered. Terra ^ 
is positive through 1200 GMT 24 April and then becomes negative 
for the remainder of the AVE 4 experiment. Its magnitude is 
consistently smaller than that of K^, except at the final time. 

To further examine variability of the two components, 
pressure-time cross sections are given in Fig. 7. A close 
correlation exists between the intense convection and the kinetic 
energy of the divergent wind component. Maximum values of 
are found in the upper troposphere at 0600 GMT 25 April when the 
second MCC is most intense. At that time, the divergent component 
comprises 5.6% of the total energy in the 300-200 mb layer. An 
upper-level maximum of also is associated with the first MCC 
at 0600 GMT 24 April. Large values of are attributable to 
large horizontal divergence at jet stream level which is hypothe- 
sized to be induced by the storm outbreak (Maddox et al. , 1981). 
Weaker, secondary maxima appear in the lower troposphere near 
the times of peak convection. Time variability of (Fig. 7) 
appears to be more closely related to the larger-scale circula- 
tion than to the storm outbreaks. Maximum values are persistently 
found at the jet stream level (300-200 mb layer). 

Contributions of the two components to total generation are 
shown in Table 7. The divergent component provides a source of 
kinetic energy to the region at every time. The maximum value 
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Pressure-time cross sections of area averaged kinetic 

A 

energy budget terms. Units of and K-, are 10 
-2 5 -o U K 

J m /100 mb and 10 J m /100 mb, respectively; 

-2 

while other units are W m /100 mb. 






P(mb) 







Table 7 . Vertical totals of the area averaged components of 
generation and horizontal flux divergence for ^ 
individual observation times. All units are W m . 


Date/Time 

->• ->• 
- v T - V * 

-v v * 


^•k^ T 

V * kV R 

VkV D 

24/00 

0.77 

- 3.19 

~ Ts 795 

- 6.07 

- 8.80 

2.73 

24/06 

- 10.48 

- 13.49 

3.01 

8.68 

2.18 

6.50 

24/12 

- 24.26 

- 26.81 

2.55 

- 1.53 

- 5.53 

4.00 

24/15 

- 12.77 

- 17.52 

4.75 

- 2.43 

- 5.75 

3.32 

24/18 

9.87 

4.94 

4.93 

- 3.35 

- 7.43 

4.08 

24/21 

5.61 

3.01 

2.60 

- 2.17 

- 6.33 

4.16 

25/00 

- 5.84 

- 7.16 

1.32 

1.03 

- 3.44 

4.47 

25/06 

2.72 

- 3.99 

6.71 

4.68 

0.46 

4.22 

25/12 

- 1.50 

- 4.90 

3.40 

15.85 

9.71 

6.14 

Composite 

< 

O 

- 7.7 

3.7 

1.6 

- 2.8 

4.4 



_o 

of 6.71 W m occurs near peak intensity of the second outbreak. 
On the other hand, the rotational wind destroys kinetic energy 
at most times. The largest negative values occur between 0600 
GMT and 1500 GMT 24 April. It is noteworthy that the two com- 
ponents generally are equally important in explaining total 
generation. This happens in spite of the fact that is much 
smaller than (Table 6). 

The pressure-time cross section for -V D * V<f> (Fig. 7) indi- 
cates only weak generation and destruction of energy until the 
second outbreak of convection. At 0600 GMT 25 April, strong 
generation of kinetic energy occurs within the 300-200 mb layer, 

and continues through the end of the experiment. Peak values of 

_o 

3.75 W m /100 mb within this layer are similar to those of the 

_o 

total generation (3.53 W m /100 mb). Figure 7 also contains 
the pressure-time cross section of Strong destruction 

of kinetic energy is evident above 400 mb throughout the AVE 4 
experiment. Maximum destruction occurs in the 100-200 mb layer 
approximately 6 h after the first MCC reaches peak intensity. 

As noted in Table 7, total generation is primarily influenced by 
the rotational wind until the second storm outbreak. 

Time variability of total horizontal flux divergence and 
its two components is given in Table 7 . The divergent wind 
produces a consistent export throughout the period. Maximum 
values of 6.50 and 6.14 W m ^/100 mb occur at 0600 GMT 24 April 
and 1200 GMT 25 April, respectively. At peak intensity of the 
first MCC (0600 GMT 24 April) , nearly 75% of the total hori- 
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zontal flux is due to the divergent wind. Over 90% of 
^■kV T is attributable to the divergent wind at 0600 GMT 25 
April when the second MCC is near peak intensity. The rota- 
tional wind, on the other hand, generally provides a source of 
energy to the region. Flux convergence is prominent during the 
period between the two convective outbreaks and at the start of 
the experiment. Relatively weak flux divergence occurs at both 
peak convection times, while a strong sink is evident at the 
final time. Even though 36 h composite values (Table 4) indi- 
cated that V-kV^ dominates the total flux, an examination of 
the individual times (Table 7) reveals that this dominance 
actually occurs only during the intense convective periods. 
Again, one should recall that is much smaller than its 
nondivergent counterpart (Table 6). 

*)■ 

Figure 7 shows pressure-time cross sections of V-kV D and 
V*kV R . Energy export due to is maximized within the 
300-200 mb layer. Greatest values at this height occur at 0600 

GMT 24 April (3.26 W m" 2 /100 mb) and at 1200 GMT 25 April (3.96 

-2 

W m /100 mb). The diagram for rotational flux shows a less 
distinct pattern with greatest values near 300 mb at the last 
observation time. 

The relative importance of the divergent wind during these 
convective outbreaks, which does not include a major cyclone, 
is similar to that noted by Chen et al. (1978) during a period 
of major cyclogensis. Important features of the two investi- 
gations include: 
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1) Term is nearly two orders of magnitude smaller than 
K d in both cases. Global studies by Chen and Wiin- 
Nielson (1976) also found small values of K^/Kp, and 
as observed here, maximum values of were found near 
jet stream level. 

2) Term V-kV^ accounts for a significant percentage of 

the total horizontal flux. In the cyclone study, 25 - 
50% of was due to the divergent wind. In the 

present case, as much as 90% of the total horizontal 
flux is attributable to during times of peak con- 
vection. 

3) In both studies, rotational generation generally is 
negative and opposes positive contributions by the 
divergent component. During the cyclone case, diver- 
gent generation was a maximum when the cyclone was 
most intense. During AVE 4, this maximum is not 
evident during the first convective outbreak, but is 
observed during the second MCC. 

4) Both components produced flux convergence at all 
levels in the cyclone case, but opposing contributions 
are observed during AVE 4. This could be due to 
differing locations and strengths of the jet streams 
with respect to the computational areas. 

5) An important conclusion is that is as important 
during the severe storm outbreaks as during the extra- 
tropical cyclone case. 
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c . Spatial distribution of energy budget terms 

Roles of the divergent and rotational wind components have 
been discussed within the context of the entire computational 
region. However, averaging over such a large area does not 
permit a detailed study of relationships between the energy 
budget terms and the enclosed convection. Therefore, spatial 
maps of selected budget terms will be presented to further 
explore interactions between the storms and their surrounding 
synoptic-scale environments. Horizontal fields of the diver- 
gent and rotational components of kinetic energy content, 
generation, and horizontal flux are presented during the life 
cycle of the second MCC (0000 GMT to 1200 GMT 25 April). 

Fields describing the surface-700 mb layer are denoted as layer 
1, while maps for the upper troposphere (400-100 mb) are de- 
noted as layer 3. Corresponding maps of the total kinetic 
energy budget have been presented by Fuelberg and Scoggins 
(1978). It should again be stressed that energetics of the 
synoptic-scale storm environments are being investigated, not 
the individual convective elements. 

1) Kinetic energy 

Spatial maps of divergent kinetic energy content are given 
in Figs. 8 and 9 for the lower and upper troposphere, respec- 
tively. At 0000 GMT, weak values of are found in both 
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layers. By 0600 GMT, however, two areas of maximum can be 
seen in the lower troposphere. The western region, extending 
from Indiana to central Texas, agrees favorably with the loca- 
tions of strongest convective activity (Fig. 5). The second 
area, over the Tennessee River Valley, is found in advance of 

the convection and it contains a broad maximum with values 

4-2 

greater than 4 x 10 J m . In the upper troposphere (Fig. 9), 

values of increase greatly between 0000 GMT and 0600 GMT. A 

4 -2 

well defined maximum greater than 25 x 10 J m is located 
over Missouri, just northwest of the MCC. This region of 
strongest K R is near the greatest upper-level velocity diver- 
gence; at 200 mb values reach 14 x 10 s At 1200 GMT, as 
the convection decreases in intensity, the field of in the 
lower troposphere (Fig. 8) appears more organized and closely 
related to storm locations. Maximum values over western 
Tennessee are located very near the tallest echo tops (Fig. 5). 
In the upper troposphere (Fig. 9), K R has decreased dramati- 
cally between 0600 GMT and 1200 GMT. Maximum values now are 
located over Illinois and Indiana, northwest of the weakening 
MCC. 

Since the rotational wind comprises a great proportion of 
the total flow, fields of K R (Figs. 10 and 11) closely resemble 

K r in the lower levels 
is maximized along a band extending from the East Coast to the 
western Gulf of Mexico (Fig. 10). In the upper troposphere at 
0000 GMT (Fig. 11), large values of extend from Texas to the 


those of Kj, (not shown). At 0000 GMT, 
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eastern Great Lakes. Maximum values are located over northern 

Texas. The developing convection, is occurring just in advance 

of the maximum (Fig. 5). By 0600 GMT (Fig. 10), a new maximum 

of forms in the lower troposphere along the eastern edge of 

the now intense complex. This maximum persists and remains 

along and east of the strongest thunderstorms through 1200 GMT. 

In the upper troposphere at 0600 GMT (Fig. 11), a new maximum 

5 -2 

has developed suddenly. Values greater than 30 x 10 J m are 
located over northern Illinois, immediately north of the convec- 
tive complex. A weaker, secondary maximum continues over northern 
Texas and Oklahoma. By 1200 GMT, strongest is located in the 
upper Ohio River Valley, well behind the MCC. The formation of 
upper-level wind maxima north of MCC's appears to be a common 
feature of these organized convective systems (Fritsch and 
Maddox, 1981a). Explanations for such development will be 
sought in following sections. 

2) Generation of kinetic energy 

Figures 12 and 13 show the spatial distribution of cross - 
contour generation by the divergent wind. At 0000 GMT, values 
of V • V* are weak in both the upper and lower troposphere. By 
0600 GMT, the lower troposphere contains destruction of energy 
near the intense convection while generation is evident to the 
east of the complex. This pattern continues through 1200 GMT. 

In the upper troposphere (Fig. 13), startling changes occur 
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Fig. 12. Horizontal maps of term (G^' 

the surface-700 mb layer (LI). Vaj 

• t 7 -2 

m W m . 
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between 0000-0600 GMT. Very strong generation of kinetic energy 
develops immediately northwest of the severe thunderstorms. 

Weaker destruction is evident south and east of the MCC. The 

_2 

intensity of the generation (greater than 120 W m ) can be 
attributed to thunderstorm induced outflow being directed pole- 

ward across the height contours toward lower values. Magnitudes 

-+ -1 
of at 200 mb are as great as 14 m s in this region. One 

should note that the area of greatest generation is located near 
the region of strong increases (Fig. 9). By 1200 GMT, strong 
generation of energy continues northwest of the complex; however, 
the intensity has decreased significantly. 

Generation by the rotational wind in the lower troposphere 
is shown in Fig. 14. Values are weak at 0000 GMT, but by 0600 
GMT an area of strong generation develops over western Tennessee 
and remains nearly stationary through 1200 GMT. Destruction of 
energy over the southern Appalachians at 0600 GMT decreases in 
intensity by 1200 GMT. In the upper troposphere (Fig. 15), 
regions of strong generation and destruction cover much of the 
AVE 4 region. At 0000 GMT, generation by V^ clearly dominates 
that due to V^ (Figs. 13 and 15); however, by 0600 GMT, genera- 
tion by the two wind components is of equal importance. An 
especially interesting feature at 0600 GMT is a couplet of 
strong generation and destruction over the upper Mississippi and 
Ohio River Valleys. Upper-level anticyclonic flow in the vicinity 
of the storms is the apparent cause for this pattern. Fritsch 
and Maddox (1981a) have noted anticyclonic flow near jet stream 
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levels of numerous MCC's. Strong generation by both components 

-2 

produces a total value greater than 200 W m over Illinois, 
near the region of strong increases (Fig. 9). Generation by 
^ persists over the upper Mississippi River Valley through 1200 
GMT; however, the area of destruction to the east decreases in 
intensity (Fig. 15). 

3) Horizontal flux of kinetic energy 

In the lower troposphere, flux convergence of energy due to 
\Tp occurs along and east of the MCC throughout its life cycle 
(Fig. 16). At peak storm intensity (0600 GMT) maximum flux 
convergence is located over Arkansas, southeast of the strongest 
thunderstorms (Fig. 5). This area moves northeastward by 1200 
GMT as does the strongest convection. Flux divergence is found 
to the rear of the convective complex at all three times. The 
upper troposphere is characterized by weak values of ^'kV^ at 
0000 GMT (Fig. 17). Flux convergence occurs over northern 
Indiana while weak flux divergence is found over the developing 
storm complex. Spectacular increases occur by 0600 GMT when 

very strong flux divergence develops over the entire convective 

_o 

complex. The maximum, containing values greater than 140 W m 
over southern Illinois, is near the region of greatest upper- 
level velocity divergence (not shown) . Strong flux convergence 
of energy has developed over northern Illinois, downwind of the 
developing Kj, maximum. One should recall that approximates 


62 





K^,; therefore, the pattern is similar to that of Fig. 11. At 
1200 GMT (Fig. 17), flux divergence continues over the weakening 
convection although values have decreased. The flux convergence 
over Illinois also has weakened. 

-y. y 

Horizontal maps of v*kV_ are given in Figs. 18 and 19. 

R 

These fields reflect the advection of Kj, patterns since the 
kV'V^ component of the flux term is near zero. Flux convergence 
generally is located east of the MCC while flux divergence 
occurs behind the storms. Greatest magnitudes are observed at 
the last observation time. In the upper troposphere at 0000 GMT 
(Fig. 19), values of V-kV are weak. By 0600 GMT, however, a 
strong center of flux divergence is located over the upper 
Mississippi River Valley, north of the convective complex and 
upwind of the K T maximum that is approximated by Fig. 11. Peak 
values of 280 W m are nearly twice those of V»kV D (Fig. 17); 
however, the two maxima are not co-located and the area due to 
is considerably smaller. Weaker flux convergence is located 
over northern Illinois and eastward. At 1200 GMT weak flux 
convergence is occurring over the convective complex that 
stretches from the Appalachians to Arkansas (Fig. 5). Maximum 
flux divergence is located over the Great Lakes, behind the wind 
maximum that is exiting the region. 
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d. Energy budgets of convection subregions 


To quantify energetics of the near storm environment, 
energy budgets were obtained for subvolumes of the total region 
that just enclosed the convective areas. These volumes, based 
only on MDR data (Table 1), moved with the convection and were 
identical to those used by Fuelberg and Scoggins (1978) . Al- 
though the limited areas were only about one-fourth the total 
computational region, it should be emphasized that their budgets 
describe energetics of the synoptic-scale flow in which the 
convection was embedded, and not the individual storms themselves. 
Figure 20 is an example of the boundaries and associated MDR 

coding for 0600 GMT 25 April. The average area of these smaller 
12 2 

subregions was 10 m . 

Time variability of the total flow energy terms is shown in 

Table 8 for the second MCC (MCC2) . Total kinetic energy increases 

5 -2 

to a peak value of 27.3 x 10 J m at 0000 GMT, then decreases 
with time through 1200 GMT. Generation due to cross -contour 
flow provides an energy source throughout the MCC's life cycle 
with a maximum value occurring near peak storm intensity (0600 
GMT). Flux divergence and negative dissipation are responsible 
for the decreases in Kp after 0000 GMT. Although flux divergence 

of energy out of the region occurs at every time after 1800 GMT, 

-2 

the maximum value of 43.9 W m is observed during greatest 
convective activity. Transfer of energy to subgrid scales of 
motion occurs at each time. As noted by Fuelberg and Scoggins 
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Table 8. 


Vertical totals of the kinetic energy budget for 
limited volumes enclosing MCC2 at individual 
observation times. All units are W m except 
Kp which is l(r J m . 


Date/Time 



-v T - 

-^ 4- 

v-kv T 

D 

24/18 

24.7 

6.0 

• — i 

-21.5 

-22.9 

24/21 

25.8 

5.1 

19.5 

11.6 

- 2.8 

25/00 

27.3 

-12.0 

4 . 6 

13.0 

- 3.6 

25/06 

25.5 

- 6.3 

52.2 

43.9 

-14.6 

25/12 

23.8 

- 7.9 

21.1 

10.3 

-18.7 

Composite 

25.4 

- 3.0 

21.0 

11.5 

-12.5 



(1978), these major energy conversions and transports near the 
convective complex are in striking contrast to the much weaker 
processes observed for the AVE 4 experiment area as a whole 
(Table 2). 

Table 9 gives components of kinetic energy content for 
MCC2. Values of Kjj show a cyclic variation that is closely 
related to the life cycle of the convection. Magnitudes in- 
crease nearly five-fold between 0000-0600 GMT as the storms 

intensify, and then decrease between 0600-1200 GMT as the storms 

S -2 

weaken. At 0600 GMT, the value of 1.72 x 10 J m constitutes 

6.7% of the total kinetic energy content. It is noteworthy that 

terms K^. and Kr decrease between 0000-0600 GMT. However, one 

should recall that a maximum of Kr did form north of the storm 

area (Fig. 11) - outside of the limited region (Fig. 20). Thus, 

placement of the limited volume is quite important. Values of 

remain negative during the MCC's life cycle, and are 

smallest at the time of peak convection. 

Generation of kinetic energy by the total flow increases 
-2 -2 

from 4.6 W m at 0000 GMT to 52.2 W m only 6 h later at peak 

storm intensity (Table 10). Although both wind components blow 

toward lower heights and thereby generate kinetic energy at 0600 

GMT, the divergent component is the primary contributor to the 

sudden increase within the limited region. Even though is a 

small fraction of the total kinetic energy content (Table 9), 

-2 

divergent generation (35.7 W m ) accounts for 68% of 
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Table 9. Vertical totals of the components of kinetic energy 

content for limited volumes enclosing MCC2j.at individual 
observation times. Units of energy are 10 J m 


Date/Time 


*R 


R D 

V*i 

24/18 

24.7 

26.3 

0.33 

-1.9 

~T73% 

24/21 

25.8 

27.4 

0.38 

-2.0 

1.5% 

25/00 

27.3 

27.6 

0.37 

-0.7 

1.4% 

25/06 

25.5 

24.0 

1.72 

0.2 

6.7% 

25/12 

23.8 

23.7 

0.63 

-0.5 

2.6% 

Composite 

25.4 

25.8 

0.69 

-1.1 

2.7% 



Table 10. Vertical totals of the components of generation and 

horizontal flux divergence for limited volumes enclosing _2 
MCC2 at individual observation times. All units are W m 


Date/Time 

4- 

-v T * V<j> 

->■ ->■ 
-v R . V<f> 

-► -► 
-v D 'V^» 

V-kV T 


V ' kV D 

24/18 

7.4 

- 6.2 

13.6 

-21.5 

-20.3 

- 1.2 

24/21 

19.5 

18.3 

1.2 

11.6 

- 0.7 

12.3 

25/00 

4.6 

7.9 

- 3.3 

13.0 

1.1 

11.9 

25/06 

52.2 

16.5 

35.7 

43.9 

5.6 

38.3 

25/12 

21.1 

14.5 

6.6 

10.3 

- 2.5 

12.8 

Composite 

21.0 

10.2 

10.8 

11.5 

- 3.3 

14.8 


. Generation decreases after 0600 GMT as the storms 

weaken. 

Horizontal flux divergence provides a major energy sink to 
the storm region, especially when the storms are most intense 
(Table 10). The divergent wind component is the major cause for 

*4* “V 

this export. After 1800 GMT, values of V-kV^ are consistently 
positive and greater in magnitude than those of V»k^ R . As 
observed with previous terms involving the divergent component, 

a major increase in V'kV^ occurs between 0000-0600 GMT. At 

_o 

0600 GMT, the divergent wind contributes 38.3 W m or 87% of 
the total flux divergence within the limited volume. The influ- 
ence of V R is much greater to the northwest of the storm region 
(Figs. 17 and 19). 

It is informative to describe vertical variations of the 
energy terms within the subvolume at peak storm intensity (0600 
GMT 25 April) . Vertical distributions of total flow energy 
terms are given in Table 11. Maximum kinetic energy content 
occurs at jet stream level (300-200 mb layer). Values are 
considerably greater than those of the composite area-time 
average (Table 2). Generation of kinetic energy due to cross- 
contour flow is a source at nearly all levels; however, greatest 
values are found in the upper troposphere. Horizontal flux 

divergence is a major exporter of energy above 500 mb, with a 

_2 

maximum of 28.8 W m /100 mb at jet level. Below 500 mb, weak 
horizontal flux convergence occurs. The vertical flux term 
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Table 11. Average kinetic energy budget for the limited area 
enclosing the ^onvection at 0600 GMT ^5 Aprjl. All 
units are W m z except Kj, which is 10 J J m . 


Pressure 
Layer (mb) 


31^/ at 

-e- 

* 

V'kV T 

3oik/3p 

D 

200-100 

4.8 

-2.2 

4.3 

7.5 

- 4.8 

- 3.8 

300-200 

6.0 

0.2 

16.9 

28.8 

-11.6 

0.5 

400-300 

4.4 

-0.1 

14.3 

10.8 

- 6.5 

-10.1 

500-400 

3.6 

-1.4 

11.0 

1.3 

2.7 

- 8.4 

600-500 

2.6 

-0.9 

4.1 

- 0.9 

7.3 

1.4 

700-600 

1.7 

-1.0 

0.5 

- 1.2 

8.1 

5.4 

800-700 

1.0 

-0.5 

-2.6 

- 0.3 

2.5 

4.3 

900-800 

0.9 

-0.3 

1.3 

- 1.1 

1.2 

- 1.5 

SFC-900 

0.5 

-0.1 

2.4 

- 1.0 

1.1 

- 2.4 

Vertical 

Total 

25.5 

-6.3 

52.2 

43.9 

0.0 

-14.6 


i 


3u)k/8p indicates strong upward transport because values of 
synoptic-scale vertical motion near the storm complex are as 
great as -16 yb s ^ at 500 mb. Transfers of energy to subgrid 
scales provide an important energy sink at 0600 GMT. Destruc- 
tion below 800 mb is due to surface friction; negative values 
also occur in the upper troposphere. Transfer of energy from 
subgrid to resolvable scales of motion (positive dissipation) is 
observed between 800-500 mb and within the 300-200 mb layer. 
Middle tropospheric positive dissipation also has been observed 
in convective storm environments by Vincent and Schlatter (1979). 

Vertical distributions of kinetic energy content at 0600 
GMT, shown in Table 12, show the importance of ^ near jet 
level. Within the 300-200 mb layer, where velocity divergence 
reaches 14 x 10 s \ the divergent wind component provides an 
11.8% contribution to K^,. An additional maximum of Kjj, occurring 
in the 700-800 mb layer, comprises 12% of JCp and is due to the 
pronounced low-level velocity convergence in the storm environ- 
ment. The vertical distribution of closely follows that of 

Kp with a maximum value at jet level. Thus, failure to include 

•+ 

Vpj in analyses of the storm region would produce energy errors 
ranging between 2-12% within individual 100 mb layers. 

Generation of kinetic energy as a function of pressure is 
shown in Fig. 21 and Table 13. Except for the 200-100 mb layer, 
the rotational component creates kinetic energy at all levels; 
maximum generation occurs near 350 mb. Within the storm volume, 
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Table 12. Components of average kinetic energy content for 
the limited area enclosing the convection at 0600 
GMT 25 April. Units of energy are 10^ J m . 


Pressure 
Layer (mb) 


Kr 

5> 

V *v 
R D 

V*T 

200-100 

4.8 

4.5 

0.22 

0.12 

4.6% 

300-200 

6.0 

5.4 

0.71 

-0.10 

11.8% 

400-300 

4.4 

4.0 

0.32 

0.08 

7.3% 

500-400 

3.6 

3.3 

0.11 

0.16 

3.1% 

600-500 

2.6 

2.5 

0.06 

0.05 

2.3% 

700-600 

1.7 

1.7 

0.08 

-0.11 

4.7% 

800-700 

1.0 

1.2 

0.12 

-0.29 

12.0% 

900-800 

0.9 

1.0 

0.07 

-0.15 

7.8% 

SFC-900 

0.5 

0.4 

0.03 

0.03 

6.0% 

Vertical 

Total 

25.5 

24.0 

1.72 

-0.21 

6.7% 



Fig. 21. Generation of kinetic energy in the liaited volume 
enclosing the convection at 0600 GMT 25 April. 

All units are W m /10Q mb. 
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Table 13. Components of average generation and horizontal flux 

divergence for the limited area enclosing 2 the convection 
at 0600 GMT 25 April. All units are W m . 


Pressure 
Layer (mb) 

-> ->• 

-V T 'V* 

4 - 

-V R -V<J> 

-V D -V(f» 

4 4 

V ' kV T 

VkV R 


200-100 

4.3 

-1.0 

5.3 

7.5 

-1.0 

8.5 

300-200 

16.9 

4.1 

12.8 

28.8 

3.8 

25.0 

400-300 

14.3 

5.0 

9.3 

10.8 

2.2 

8.6 

500-400 

11.0 

3.5 

7.5 

1.3 

0.7 

0.6 

600-500 

4.1 

1.5 

2.6 

- 0.9 

0.3 

- 1.2 

700-600 

0.5 

0.5 

0.0 

- 1.2 

0.3 

- 1.5 

800-700 

- 2.6 

0.3 

- 2.9 

- 0.3 

0.2 

- 0.5 

900-800 

1.3 

1.5 

- 0.2 

- 1.1 

-0.5 

- 0.6 

SFC-900 

2.4 

1.1 

1.3 

- 1.0 

-0.4 

- 0.6 

Vertical 

Total 

52.2 

16.5 

35.7 

43.9 

5.6 

38.3 


the divergent component provides a much stronger energy source 
above 500 mb than does the rotational component. Within these 
layers, over 75% of is due to V^. Near 750 mb, however, 

Vp predominantly crosses the contours toward higher values, 
thereby destroying kinetic energy. Since generation by the 
divergent wind greatly exceeds that due to the nondivergent 
wind, the profile for total generation is similar to that of 

4 4* (Fig- 21). 

Vertical profiles of the horizontal flux terms at 0600 GMT 
(Fig. 22 and Table 13) show that greatest outflow from each com- 
ponent occurs near the jet stream level (250 mb); however, 

-y -*■ 

magnitudes from far exceed those from V^. Within the 300-200 
mb layer, nearly 87% of the total flux divergence is due to V^. 
Neither component produces large magnitudes in the lower and 
middle troposphere. These results indicate that failure to 
include in computational work would produce serious misrepre- 
sentations of kinetic energy generation and transport within the 
storm environment. 

To verify the role of the divergent wind during intense 
convective outbreaks, the energetics of the first MCC (MCC1) 
which occurred early in the AVE 4 experiment (0000-1500 GMT 24 
April) are presented. Boundaries for the subvolumes that enclose 
it again are identical to those of Fuelberg and Scoggins (1978). 
One should note the numerous similarities between these two 
cases. Values of for MCC1 (Table 14) show a cyclic variation 
which is similar to that observed for MCC2 (Table 9). A maximum 
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Fig. 22. Horizontal flux divergence of kinetic energy 

in the limited volume enclosing the convection 
at 0600 GMT 25 April. All units are W ■ /100 mb. 
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Table 14. Vertical totals of the components of kinetic energy 
content for limited volumes enclosing MCC1 at 
individual observation times. Units of energy are 
1(T J m . 


Date/Time 

s 

*R 

“d 

V -V 
R D 

V*T 

24/00 

30.8 

31.2 

0.42 

-0.8 

1.4% 

24/06 

24.1 

23.6 

1.00 

-0.5 

4.1% 

24/12 

26.4 

25.4 

0.87 

0.1 

3.3% 

24/15 

25.2 

25.5 

0.36 

-0.6 

1.4% 

Composite 

26.6 

26.4 

0.66 

-0.5 

2.5% 



5 -2 

value of 1 x 10 J m , or 4.1% of Kp, is observed at peak storm 
intensity (0600 GMT 24 April) . The value of doubles between 
0000-0600 GMT, then decreases after 0600 GMT as the storms 
weaken. As with MCC2, values of Kp and K R decrease within the 
subvolume as the storms intensify; however, maxima of K R and Kp 
again form north and outside of the storm volume. 

Generation of kinetic energy by the total flow increases 
from 30.2 W m~ 2 at 0000 GMT to 47.5 W m 2 6 h later, at peak 
storm intensity (Table 15). Contrary to results for MCC2 (Table 
10) , the increase is due mostly to variations in rotational 
generation. Actually, however, horizontal fields of in 

the 400-100 mb layer (not shown) indicate strong generation of 
kinetic energy to the north of MCC1 and destruction to the 

south. This pattern is similar to that observed with MCC2 (Fig. 

-2 

13). Maximum values of 104 W m also are comparable to those 
of MCC2 (Fig. 13). With MCC1, however, the greatest effect of ^ 
occurs outside the limited region, thus accounting for the com- 
paritively small values of divergent generation in Table 15. 

Thus, some of the contrasts between the two systems are attri- 
butable to the relative locations of the subvolumes with respect 
to particular energy features. Both total and rotational genera- 
tion decrease dramatically between 0600-1200 GMT as the storms 
weaken (Table 15). Spatial maps of divergent generation show 
that it too decreases after 0600 GMT. 

As with MCC2, horizontal flux divergence provides a major 
energy sink during the life cycle of the first convective com- 
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Table 15. Vertical totals of the components of generation and 
horizontal flux divergence for limited volumes 
enclosing MCC1 at_ jndividual observation times. 

All units are W m 



Date/Time 

-v T ‘ V* 


-v D - V 4 » 

-4- ->• 

V'kV T 

’• k \ 

^. k ? D 


24/00 

30.2 

15.5 

14.7 

- 1.2 

0.3 

- 1.5 

24/06 

47.5 

40.2 

7.3 

44.7 

19.7 

25.0 

24/12 

11.0 

2.2 

8.8 

28.4 

7.8 

20.6 

24/15 

17.3 

9.9 

7.4 

8.7 

- 1.7 

10.4 

Composite 

26.5 

17.0 

9.5 

20.1 

6.5 

13.6 


plex (Table 15). Again, a significant increase in 

occurs between 0000-0600 GMT 24 April as the complex develops. 

The nondivergent component also increases during this 6 h period, 
and both components contribute nearly equally to ?»kV T during 
peak convective activity. 

Vertical distributions of kinetic energy content at 0600 
GMT 24 April (Table 16) show that both components exhibit maxima 
at jet stream level. Within the 300-200 mb layer, for MCC1 
is only one-third that of MCC2 (Table 12) , although it still 
contributes 4.4% of K,p. Spatial maps of for the upper tropo- 
sphere (not shown) indicate maximum values occurring on the 
poleward side of the convection, similar to that found for MCC2 
(Fig. 9). A secondary maximum is observed in the lower tropo- 
sphere near 850 mb (Table 16). 

Table 17 gives vertical distributions for the components of 
generation and horizontal flux at 0600 GMT 24 April. Maximum 
generation occurs in the upper troposphere above 300 mb and 
appears due to the nondivergent wind. Again, it should be noted 
that the generation maximum due to occurs north of the sub- 
volume and therefore does not prominently appear in the Table. 
Both wind components produce horizontal flux divergence of 
energy above 500 mb and flux convergence below 500 mb. Hori- 
zontal flux due to Vjj dominates V*kV^, at nearly every level 
during peak convective activity. 
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Table 16. Components of average kinetic energy content for 
the limited area enclosing the convection a£ 0600 
GMT 24 April. Units of energy are 10 J m . 


Pressure 
Layer (mb) 

*T 

Kr 

5> 

V • V 
R D 

VS 

200-100 

5.2 

5.1 

0.10 

-0.04 

1.9% 

300-200 

5.4 

5.5 

0.24 

-0.30 

4.4% 

400-300 

3.7 

3.7 

0.07 

-0.07 

1.9% 

500-400 

3.0 

2.6 

0.06 

0.32 

2.0% 

600-500 

2.4 

2.1 

0.07 

0.27 

2.9% 

700-600 

1.7 

1.6 

0.11 

-0.01 

6.5% 

800-700 

1.3 

1.4 

0.13 

-0.19 

10.0% 

900-800 

1.0 

1.2 

0.17 

-0.34 

17.0% 

SFC-900 

0.4 

0.4 

0.09 

-0.12 

22.5% 

Vertical 

Total 

24.1 

23.6 

1.0 

-0.5 

6.1% 



Table 17. Components of average generation and horizontal flux 

divergence for the limited area enclosing^the convection 
at 0600 GMT 24 April. All units are W m 


Pressure 
Layer (mb) 

-)■ 

-v v * 

-V R * V<t> 

-v D *v<j) 

4 

V-kV T 

■+ 4- 

V^R 

VkV D 

200-100 

12.2 

11.0 

1.2 

10.9 

4.3 

6.6 

300-200 

11.8 

10.0 

1.8 

29.5 

10.4 

19.1 

400-300 

6.2 

4.6 

1.6 

12.8 

5.2 

7.6 

500-400 

4.7 

2.1 

2.6 

3.7 

2.7 

1.0 

600-500 

3.3 

1.8 

1.5 

- 2.3 

- 0.2 

- 2.1 

700-600 

4.4 

3.1 

1.3 

- 3.5 

- 0.6 

- 2.9 

800-700 

2.9 

3.6 

-0.7 

- 3.2 

- 1.2 

- 2.0 

900-800 

0.5 

2.6 

-2.1 

- 2.3 

- 0.9 

- 1.4 

SFC-900 

1.5 

1.4 

0.1 

- 0.9 

0.0 

- 0.9 

Vertical 

Total 

47.5 

40.2 

7.3 

44.7 

19.7 

25.0 


The major point to note is that the energetics of MCC1 
generally are similar to those of MCC2. The divergent wind 
component is very important in each case. These agreements lend 
credibility to the overall study. 
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SUMMARY AND CONCLUSIONS 


Relative contributions of divergent and rotational wind 
components to the kinetic energy balance of two severe storm 
environments have been investigated. Synoptic-scale rawinsonde 
data at 3 or 6 h intervals were available from NASA's AVE 4 
experiment. Energy budgets were obtained for the entire com- 
putational region and for subvolumes that barely enclosed and 
moved with the two convective complexes . This study was an 
extension of the total kinetic energy analysis of AVE 4 by 
Fuelberg and Scoggins (1978). 

The divergent wind component was found to be very important 
in the synoptic-scale environments of the two mesoscale convec- 
tive complexes. Although usually comprised less than 10% of 

— >■ 

the total kinetic energy, was a major cause for energy gen- 
eration by cross-contour flow and for horizontal flux divergence. 
Within the subvolumes, as much as 68% of the total cross-contour 
generation and 87% of the total horizontal flux in the upper 
troposphere was due to during peak convective periods. The 
time variability of divergent energy transformations displayed a 
close correlation to the intensity of the convection. Pressure- 
time cross sections of K^, -V^‘^<J>, and ^*kV^ each 
have indicated distinct maxima at times of peak convection. 
Since strong synoptic-scale, upper-level divergence may be a 
common by-product of certain types of storm regions (Fritsch and 
Maddox, 1981a and b) , the importance of during AVE 4 may be 
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characteristic of other cases as well. A comparison of current 
results with those from a period of cyclogenesis (Chen et al. , 
1978) suggests that in major convective areas may affect the 
general circulation as much or more per unit area than in mature 
cyclone situations. 

Generation of kinetic energy by the divergent component 
appears to have been a major factor leading to the creation of 
upper-level wind maxima north of the storm areas. Recent inves- 
tigations suggest that similar wind perturbations are common 
occurrences near certain intense storm complexes (Fritsch and 
Maddox, 1981a). Thus, further investigation of the role of the 
divergent wind component may be a fruitful method for exploring 
the effects of convection on the larger-scale environment. 

As noted by Chen et al. (1978) current numerical prediction 
models generally reduce divergence (and therefore the divergent 
wind) in order to prevent spurious effects due to gravity waves. 
Current results, however, suggest that failure to adequately 
represent the divergent wind could lead to serious misrepresen- 
tations of the kinetic energy balance and possible errors in 
derived forecasts. 

The results of this study suggest several areas for future 
investigation. Current findings should be verified with those 
from other storm outbreak cases. The use of subsynoptic scale 
data sets, such as those from AVE-SESAME *79, would be especially 
useful. In addition, the conversion of kinetic energy between 
rotational and divergent components should be explored to better 
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understand the formation of wind maxima near the storms. 
Krishnamurti and Ramanathan (1982) have shown that this process 
is important in producing monsoon circulations near India. The 
energetics team at Saint Louis University is actively pursuing 
these follow up studies. 
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APPENDIX 


Sensitivity Analysis 

Sensitivity of the energy budget to errors in the input 
rawinsonde data was evaluated quantitatively. Only the effects 
of random errors were considered, i.e., systematic errors and 
computational limitations were not investigated. 

In several previous studies, random perturbations simu- 
lating rawinsonde errors have been introduced into the original 
station soundings. "Perturbed energy budgets" then were com- 
puted for comparison with original values (Vincent and Chang, 
1975; Ward and Smith, 1976; Fuelberg and Scoggins, 1980; 
Robertson and Smith, 1980; Jedlovec and Fuelberg, 1981). In a 
recent investigation, Belt and Fuelberg (1982) documented the 
ability of their objective analysis and filtering procedures to 
decrease magnitudes of deliberately introduced sounding per- 
turbations. They noted that standard deviations of height 
perturbations were reduced by approximately 30% at the grid 
points, while the more strongly filtered wind perturbations were 
reduced by approximately 50%. Since the objective analysis 
procedures of Belt and Fuelberg are almost identical to those of 
the present study (Section 3-b), an alternative technique for 
sensitivity analysis was devised that is based on their findings. 

Our procedure was to introduce computer-generated random 
perturbations at the 42 computational grid points closest to 
corresponding rawinsonde stations. The perturbations then were 
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filtered horizontally so that effects of the "errors" were 
spread to surrounding grid points. By carefully selecting 
standard deviations of the input perturbations and their degree 
of filtering, it was possible to closely simulate the way that 
the Barnes objective analysis technique would have treated 
perturbations if they had been introduced at the stations them- 
selves. Table 18 gives standard deviations for the normally 
distributed perturbations appearing in the final data grids for 
the current sensitivity study. These grid point values generally 
are either 30% or 50% smaller than station values proposed by 
Kurihara (1961) because current values reflect the ability of 
computational procedures to filter and thereby partially remove 
sounding errors. 

The 0600 GMT 25 April observation was chosen for study 
since it was the time of maximum convective activity and greatest 
upper-level divergence. Ten sets of the 0600 GMT data, each 
with a different combination of perturbations at the individual 
50 mb data levels, were made. The ten different sets should 
represent a good sample of situations that can be encountered 
with actual data. Figure 23 is an example of height perturba- 
tions at 200 mb. One should note that individual grid point 
deviations are as large as 75 m, considerably greater than 
values given in Table 18. This is possible because a normal 
distribution was allowed. Standard deviations (Table 18) and 
visual appearances of the height (Fig. 23) and wind perturbation 
(not shown) patterns show good agreement with those of Belt and 
Fuelberg (1982) and Jedlovec and Fuelberg (1981). Thus, 
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Table 18. Standard deviations of normally distributed grid 
point perturbations 


Pressure 

Level 

(mb) 

Wind 

Direction 

(deg) 

Wind 
Speed, 
(m s L ) 

Height 

(m) 

100 

5.3 

2.5 

30.2 

200 

4.5 

2.2 

24.2 

300 

3.4 

1.8 

15.5 

500 

6.1 

1.0 

8.1 

700 

1.5 

0.7 

7.0 


900 


9.2 


0.7 


3.1 




our sensitivity procedure should be a valid simple 
alternative to those cited initially. Since objective analysis 
is not required, the current procedure saves considerable computer 
time. Once the perturbed data grids were prepared, all subsequent 
procedures for computing the divergent wind, vertical motion, 
energy budgets, etc., were identical to those described in 
Section 3-b. 

Mean absolute differences between the original area aver- 
aged budget at 0600 GMT and the ten budgets derived from the 
perturbed data are given in Table 19 for the entire region and 
in Table 20 for the limited volume enclosing MCC2 (Fig. 20). 
Results indicate that term sensitivity increases with altitude 
because the rawinsonde data were assumed to become less accurate 
at the higher levels. At any particular level, one can be more 
confident of the signs of terms having large values than those 
near zero. As noted by Jedlovec and Fuelberg (1981), mean 
absolute differences for the smaller area (Table 20) are some- 
what larger than those for the entire region (Table 19); however, 
in most cases original values for the smaller area also are 
greater . 

Quantities involving the divergent wind were somewhat less 
sensitive to data errors than expected. Belt and Fuelberg 
(1982) observed that strong, well defined patterns of horizontal 
divergence were not greatly affected by data errors. On the 
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Table 19. Area averaged kinetic energy budget and mean absolute differences 

(parentheses) for the entire computational region at 0600 GMT 25 April. 
Values of egergy content are in 10 J J m while others 
are in W m . 


Pressure 

Layer 

(mb) 


*R 

*D 

V • V 
R D 

aKj/at 

3o)k/3p 

D 

400-100 

13.10 

(0.17) 

12.80 

(0.14) 

0.56 

(0.04) 

-0.20 

(0.02) 

0.4 

(0.2) 

-3.5 

(0.3) 

-2.8 

(1.3) 

700-400 

4.30 

(0.05) 

4.12 

(0.05) 

0.15 

(0.01) 

0.02 

(0.0) 

-2.4 

(0.2) 

2.4 

(0.3) 

3.0 

(0.5) 

SFC-700 

1.45 

(0.02) 

1.41 

(0.01) 

0.14 

(0.0) 

-0.09 

(0.0) 

-0.5 

(0.0) 

1.1 

(o.o) 

-0.8 

(0.1) 

SFC-100 

18.90 

(0.18) 

18.30 

(0.16) 

0.86 

(0.05) 

-0.26 

(0.02) 

-2.6 

(0.2) 

0.0 

(0.0) 

-0.6 

(1.4) 


-V T ‘V(J) 

-y -y 

-v R *v«t> 


+ + 
V-kV T 




400-100 

2.7 

(1.2) 

- 2.9 
(0.7) 

5.6 

(0.5) 

3.1 

(0.7) 

-1.1 

(0.5) 

4.1 

(0.3) 


700-400 

- 1.8 
(0.3) 

- 2.4 
(0.2) 

0.7 

(0.2) 

1.2 

(0.1) 

1.0 

(0.1) 

0.2 

(0.0) 


SFC-700 

1.8 

(0.1) 

1.4 

(0.1) 

0.4 

(0.1) 

0.4 

(0.0) 

0.5 

(0.0) 

-0.1 

(0.0) 


SFC-100 

2.7 

(1.2) 

- 4.0 
(0.8) 

6.7 

(0.6) 

4.7 

(0.7) 

0.5 

(0.5) 

4.3 

(0.3) 



Table 20. 


Area averaged kinetic energy budget and mean absolute differences 
(parentheses) for the limited a rea,. enclosing MCC2 at 0600 GMT 25 April. 
Values of energy content are in lO'' J m while others 
are in W m 


Pressure 


Layer 

(mb) 

*T 

St 

*D 

4 “4 

V • V 
R D 

aK^/at 

8o)k/3p 

D 

400-100 

15.20 

(0.32) 

13.80 

(0.26) 

1.25 

(0.06) 

0.11 

(0.08) 

-2.1 

(0.6) 

-22.8 

(1.2) 

-13.5 

(4.0) 

700-400 

7.92 

(0.13) 

7.57 

(0.13) 

0.26 

(0.01) 

0.10 

(0.03) 

-3.2 

(0.6) 

18.0 

(1.1) 

- 1.5 
(2.2) 

SFC-700 

2.41 

(0.04) 

2.61 

(0.03) 

0.22 

(0.01) 

-0.41 

(0.01) 

-0.8 

(0.1) 

4.8 

(0.2) 

0.6 

(0.2) 

SFC-100 

25.50 

(0.36) 

24.00 

(0.34) 

1.73 

(0.07) 

-0.21 

(0.10) 

-6.2 

(0.7) 

0.0 

(0.0) 

-14.4 

(5.2) 


-)• 4 - 

-V T -V<|> 

\-h 

■+ -4 

-v d -v4> 

4 - •> 

V'kV T 


*' k *D 


400-100 

35.5 

(2.9) 

8.0 

(1.8) 

27.5 

(1.4) 

47.0 

(4.5) 

5.0 

(2.1) 

42.0 

(2.8) 


700-400 

15.6 

(0.7) 

5.4 

(0.3) 

10.2 

(0.5) 

- 0.8 
(0.4) 

1.3 

(0.3) 

- 2.1 
(0.3) 


SFC-700 

1.1 

(0.2) 

2.9 

(0.1) 

- 1.8 
(0.2) 

- 2.3 
(0.2) 

-0.7 

(0.1) 

- 1.7 
(0.1) 


SFC-100 

52.2 

(3.3) 

16.3 

(1.8) 

35.9 

(1.7) 

43.9 

(4.3) 

5.7 

(2.0) 

38.2 

(2.9) 




other hand, weaker and less organized patterns were altered more 

->■ 

significantly. Thus, in the present study, terms involving 
would be smaller in magnitude and more sensitive to error at 
times other than 0600 GMT 25 April when 200 mb divergence reached 
14 x 10" 5 s' 1 . 

The sensitivity of field patterns was evaluated quanti- 
tatively by computing linear correlation coefficients between 
the original and ten perturbed budgets. Values were almost 
always greater than 0.95 for either area. Horizontal fields of 
several terms involving the divergent wind are given in Fig. 24 
for the 400-100 mb layer. For each term, the diagram depicts 
the pattern having the lowest correlation of the ten runs. 
Corresponding original fields are shown in Figs. 9, 13 and 17. 
Generally, patterns for these "worst" cases show excellent 
agreements with their originals; however, magnitudes of some 
individual grid point values occasionally are altered signifi- 
cantly. For example, maximum grid point deviations in the 
400-100 mb layer were 7.1 x 10^ J m 2 for K^, 49.3 W m ^ for 
HF d , and 21.2 W m" 2 for G D< 

In summary, random errors in rawinsonde data generally 
should not affect the interpretations of this study. Good space 
and time continuities of spatial patterns provide additional 
confidence . 
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convection within the area. Although Kp usually comprised less than 10 percent of 
the total kinetic energy content within the storm environment, as much as 87 percent 
of the total horizontal flux divergence and 68 percent of the total cross-contour 
generation was due to the divergent component in the upper atmosphere. Generation 
of kinetic energy by the divergent component appears to have been a major factor in 
the creation of an upper- level wind maximum on the poleward side of one of the com- 
plexes. A random error analysis is presented to assess confidence limits in the 
various energy parameters. 
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